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REACTION  OF  MICROPARTICLES  BY  THE  DIFFUSION 
OF  REACTIVE  GASES  THROUGH  POROUS  SHELLS 


1 .  INTRODUCTION 

Recently,  we  developed  a  new  technique  to  study  the 
heterogeneous  reaction  between  acid  aerosols  and  basic  gases.1 
The  technique,  single  particle  electrodynamic  balance  (SPEB) , 
uses  both  oscillating  and  static  electric  fields  to  stabilize 
single,  charged  droplets  at  a  fixed  point  in  space.  The  droplet 
mass  is  determined  from  the  weight-balancing  static  electric 
field  that  maintains  the  droplet  at  the  focal  plane  of  a 
telemicroscope.  By  monitoring  the  time-dependent,  weight¬ 
balancing  electric  field,  we  can  analyze  processes  that  cause 
changes  in  the  droplet  mass.  Rubel  and  Gentry1  studied  the 
heterogeneous  reaction  between  single  phosphoric  acid  droplets 
and  ammonia  gas  by  recording  the  increasing  weight-balancing 
voltage  during  droplet  reaction.  We  found  that  the  droplet 
reactions  are  sequentially  controlled  by  surface  phase,  gas  phase 
diffusion-controlled,  and  porous  shell  diffusion-controlled 
reactions.  Observing  the  droplet  through  a  telemicroscope,  we 
noted  that  gas  phase  diffusion-controlled  reactions  terminate 
after  particle  crystallization. 

The  time  of  particle  crystallization  decreased  with 
increasing  ammonia  gas  partial  pressure  and  decreasing  droplet 
size.  This  correlation  supported  the  hypothesis  that  the 
development  of  a  critical  supersaturation  of  ammonium  phosphate 
initiates  phosphate  crystallization  in  the  droplet.  To  test  the 
hypothesis,  we  derived  a  model  for  the  time-dependent  ammonium 
phosphate  saturation  field  inside  the  droplet.*  Calculations 
showed  that,  for  a  wide  range  of  experimental  conditions,  all 
droplets  were  characterized  by  the  same  surface  supersaturation 
of  ammonium  phosphate  at  the  time  of  crystallization.  This 
result  substantiated  the  hypothesis  that  a  critical  supersatur¬ 
ation  of  ammonium  phosphate  initiates  the  formation  of  a  porous 
shell  at  the  droplet  surface.  Further  evidence  for  the  formation 
of  the  surface  shells  was  found  in  photographic  records  that 


1Rubel,  G.O.,  and  Gentry,  J.W. ,  "Investigation  of  the  Reaction 
Between  Single  Aerosol  Acid  Droplets  and  Ammonia  Gas," 

J.  Aerosol  Science  Vol.  15,  No.  6,  p  661  (1984). 

*Rubel,  G.O.,  and  Gentry,  J.W.,  Unpublished  data,  1987. 


showed  the  appearance  of  opaque  annuli  at  the  droplet  surface  at 
the  time  of  crystallization. 

The  formation  of  the  porous  shell  terminates  the  gas 
phase  diffusion-controlled  reactions,  and  porous  shell  diffusion- 
controlled  reactions  became  rate  limiting.  This  study  was 
conducted  to  develop  a  semiempirical  model  that  describes  the 
porous  shell  diffusion-controlled  reaction  of  layered  micro¬ 
particles.  In  particular,  we  derive  a  model  that  predicts  the 
porous  shell  diffusion-controlled  reaction  for  phosphoric  acid 
droplets  with  surface  layers  of  ammonium  phosphate.  The  porous 
shell  reaction  model  is  parameterized  with  an  effective  diffusion 
coefficient  to  describe  the  transport  of  ammonia  gas  through  the 
ammonium  phosphate  shell.  Further  depletion  of  the  phosphoric 
acid  reactant  inside  the  droplet  is  modeled  after  a  shrinking- 
core  reaction  model.  The  porous  shell  reaction  model  is  compared 
to  data  obtained  from  the  SPEB.  We  determined  the  characteristic 
pore  radius  of  the  porous  shell  as  a  function  of  reaction 
conditions  using  the  parallel  pore  model  of  Wheeler.2 

2.  REACTION  OF  ACID  DROPLETS  IMMERSED  IN  AMMONIA  GAS 

The  use  of  the  SPEB  technique  to  measure  the  hetero¬ 
geneous  reaction  dynamics  of  single  acid  droplets  has  been 
described  in  detail  earlier.1  The  central  component  of  the 
experimental  system  is  the  electrodynamic  balance  that  stabilizes 
single,  charged  droplets  at  the  null  point  of  an  ac-driven 
hyperboloidal  electric  field.  The  droplet  weight  is  balanced 
against  a  static  electric  field  that  maintains  the  droplet  at  the 
focal  plane  of  a  telemicroscope.  The  heterogeneous  reaction  rate 
of  the  single  droplet  is  determined  from  the  rate  at  which  the 
balancing  electric  field  increases  with  time.  Using  a  35mm 
objective,  the  system  has  the  additional  capability  to  monitor 
the  time-dependent  morphology  of  the  reacting  droplet.  In  the 
earlier  study,  Rubel  and  Gentry1  measured  the  heterogeneous 
reaction  dynamics  of  phosphoric  acid  droplets  immersed  in  ammonia 
gas.  The  investigation  involved  droplet  sizes  between  42  and 
72  n m  and  ammonia  gas  partial  pressures  between  115  and 
1000  dyn/cm2.  All  experiments  were  conducted  at  one  atmosphere 
of  pressure. 

Figure  1  shows  the  extent  of  reaction  of  phosphoric  acid 
droplets  measured  as  a  function  of  time  in  the  SPEB.  The  extent 

2Wheeler,  A.,  Catalysis .  Vol.  2,  p  105,  Reinhold,  New  York,  1955. 


Figure  1.  Comparison  of  Measured  Extent  of  Reaction  (0)  and 

Theoretical  Surface  Phase  (S)  and  Gas  Phase  Diffusion- 
Controlled  (D)  Extent  of  Reaction  for  Phosphoric  Acid 
Droplets  Immersed  in  Ammonia  Gas. 
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of  reaction  is  the  ratio  of  the  number  of  ammonia  molecules 
reacted  with  the  droplet  to  the  number  of  phosphoric  acid 
molecules  initially  in  the  droplet.  This  follows  from  the  1:1 
stoichiometry  of  the  reaction.  Four  separate  reactions  are 
shown,  each  characterized  by  different  ammonia  gas  partial 
pressures.  The  data  (0)  is  compared  to  theoretical  surface  phase 
(S)  and  gas  phase  diffusion-controlled  (D)  reaction  rates. 1 
Thus,  the  reaction  dynamics  are  initially  controlled  by  surface 
phase  kinetic  reactions  and  then  by  gas  phase  diffusion- 
controlled  reactions.  Gas  phase  diffusion-controlled  reactions 
terminate  after  particle  crystallization.  The  onset  of  particle 
crystallization  is  monitored  with  a  telemicroscope. 

3.  ONSET  OF  PARTICLE  CRYSTALLIZATION 

Figure  1  shows  that  the  time  of  particle  crystallization 
decreased  with  increasing  ammonia  gas  partial  pressure  or 
equivalently  increasing  reaction  rate.  For  the  runs  shown  in 
Figure  1,  the  time  of  particle  crystallization  decreased  from 
20  to  3.5  sec  as  the  ammonia  gas  partial  pressure  increased  from 
115  to  1000  dyn/cm2.  In  a  later  paper  (unpublished  data,  1987), 
we  derived  a  quantitative  model  that  predicts  the  ammonium 
phosphate  saturation  field  inside  the  reacting  droplet.  The 
experimenters  assumed  that  the  reaction  product,  ammonium 
phosphate,  obeyed  Fick's  law  of  diffusion  and  was  generated  at 
the  droplet  surface  at  a  rate  given  by  gas  phase  diffusion- 
controlled  transport  theory.  We  presented  a  closed  form  solution 
for  the  internal  droplet  phosphate  saturation  field  as  a  function 
of  time  and  ammonia  gas  partial  pressure.  Figure  2  shows  the 
phosphate  saturation  fields  computed  at  particle  crystallization 
times  of  20,  7,  4.5,  and  3.5  sec,  corresponding  to  the  reactions 
shown  i.  Figure  1.  At  particle  crystallization,  the  surface 
supersaturation  of  ammonium  phosphate  is  1.88,  1.94,  2.02,  and 
1.90  for  the  ammonia  gas  partial  pressures  of  115,  400,  800,  and 
1000  dyn/cm2,  respectively.  Rubel  and  Gentry  conclude  that  the 
development  of  a  critical  supersaturation  of  ammonium  phosphate 
started  particle  crystallization. 

4 .  THEORY 

To  develop  a  porous  shell  diffusion-controlled  reaction 
model,  we  propose  the  following  mechanism  to  describe  the  state 
of  the  droplet  during  reaction.  At  the  time  of  crystallization, 
all  droplet  surfaces  are  characterized  by  the  same  value  of  the 
ammonium  phosphate  supersaturation.  The  excess  free  energy 
represented  by  the  phosphate  supersaturation  provides  the  driving 
force  for  the  phase  transition.  The  phase  transition  relieves 
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the  phosphate  supersaturation,  resulting  in  the  formation  of  a 
porous  shell  at  the  droplet  surface.  In  addition,  an  even  larger 
annulus,  of  saturated  ammonium  phosphate  is  circumscribed  by  the 
porous  shell.  Further  diffusion  of  ammonia  gas  through  the  shell 
results  in  the  supersaturation  of  this  phosphate-saturated 
annulus  and  the  continuous  build-up  of  the  porous  shell.  Thus, 
after  crystallization,  the  droplet  reactions  are  controlled  by 
the  diffusion  of  ammonia  gas  through  the  porous  shell  that  grows 
as  the  acid  core  shrinks. 

To  model  the  porous  shell  diffusion-controlled  reaction 
dynamics,  the  salient  properties  of  the  porous  shell  must  be 
specified,  that  is,  its  time-dependent  thickness  and  the  ammonia 
gas  concentration  field  within  the  porous  shell. 


REDUCED  RADIAL  COORDINATE  r/R 


Figure  2.  Theoretical  Ammonium  Phosphate  Saturation  Field  Inside 
a  Phosphoric  Acid  Droplet  Reacting  with  Ammonia  Gas 
(Ammonia  gas  pressure  P  is  in  units  of  dyn/cm2) . 


4 . 1  Shell  Thickness . 

At  crystallization,  all  ammonium  phosphate  exceeding  the 
saturation  concentration  cs  precipitates  from  the  solution  into 
the  body  of  the  porous  shell.  As  shown  in  Figure  2,  the  super¬ 
saturation  zone  extends  inward  from  the  droplet  surface  about 
one-tenth  the  droplet  radius.  The  thickness  through  which  the 
concentration  field  exceeds  saturation  is  R  -  8  ,  where  5  is  the 
thickness  of  the  supersaturated  zone  (Figure  3a) .  The  mass  of 
phosphate  that  precipitates  as  a  porous  shell  at  crystallization 

mph  is 

R 

mph  =  /tc<r)  “  cs]4,rr2dr  (1) 

8 

where  c(r)  is  the  radially  dependent  phosphate  concentration  at 
crystallization.  The  thickness  of  the  shell  at  crystallization 
r c  is  defined  by  the  relation 

4  ir  R3/3  -  4t(R  ~tc)3/3  =  mph/ ( 1  ~  «  )/>  ph  (2) 

where  />ph  is  the  bulk  density  of  ammonium  phosphate,  and  t  is  the 
porosity  of  the  porous  shell.  In  these  studies,  the  shell 
thickness  is  much  less  than  the  droplet  radius.  Substituting 
equation  1  into  equation  2,  the  shell  thickness  at  crystal¬ 
lization  can  be  written  as 

R 

rc  =  R  -  {r3  -  3cs  J[S(r)  -  l]r2/[pph(l  -  «)]  dr}1/3  (3) 

8 


Figure  3.  Schematic  of  Ammonium  Phosphate  Saturation  Field 
before  Crystallization  (a)  and  Ammonia  Gas 
Concentration  Field  Inside  c^  and  Outside  cex  the 
Porous  Shell  (b) . 


where  S(r)  is  the  phosphate  saturation  ratio.  The  saturation 
ratio  is  defined  as  S(r)=  c(r)/cs,  cs  being  the  phosphate 
saturation  concentration  in  the  acid  solution.  As  the  reaction 
proceeds,  the  shell  thickens.  The  thickness  of  the  shell  is 
related  to  the  extent  of  droplet  reaction  as  follows.  The  number 
of  phosphate  molecules  np^  in  the  droplet  at  time  t  is  expressed 
as 

nph(t)  =  nph(tc)  +  47TCp{(R  -tc)3  -  [R  -r(t)]3}/3  (4) 

where  n  h(t)  is  the  sum  of  the  phosphate  molecules  in  the  droplet 
at  crystallization  riph(tc)  and  the  phosphate  molecules  created  by 
the  reaction  between  phosphoric  acid  molecules  and  ammonia  gas 
diffusing  through  the  porous  shell.  Here  cps  is  the  molar  volume 
concentration  of  phosphoric  acid  molecules  in  the  saturated 
annulus.  We  define  the  extent  of  reaction  &  as  &  s  nph/np(°)# 
where  np(0)  is  the  initial  number  of  moles  of  phosphoric  acid  in 
the  droplet.  Then,  the  time-dependent  shell  thickness  is  related 
to  the  extent  of  reaction  by 

r(t)  =  R  -  {(R  -Tc)3  -  3np  ( 0)  [  &  (t)  -  £c]/(4irCp)J  1/3  (5) 

4 . 2  Ammonia  Concentration  Field. 

The  ammonia  gas  concentration  inside  the  porous  shell  is 
shown  qualitatively  in  Figure  3b.  The  ammonia  gas  concentration 
at  the  liquid  acid  core  surface  is  set  to  zero  because  the 
reaction  between  the  acid  and  ammonia  molecules  is  rapid  compared 
to  the  porous  shell  diffusion  rate.  At  the  droplet  surface,  a 
discontinuity  exists  in  the  slope  of  the  concentration  profile 
owing  to  the  different  boundary  conditions  imposed  onto  the 
interior  c^  and  exterior  cex  solutions  for  the  ammonia  gas 
concentration  field.  We  establish  the  ammonia  gas  concentration 
within  the  porous  shell  by  forming  an  imaginary  sphere  of  radius 
r  so  that  its  surface  lies  within  the  porous  shell.  The 
transport  rate  of  ammonia  gas  through  the  shell  is 

47rr2De8ci/3r 
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where  De  is  the  effective  porous  shell  diffusion  coefficient. 
Under  steady  state  conditions,  the  flux  is  independent  of  radius 
and 

4nr2Deac^/3r  =  N  =  constant  (6) 


Integrating  equation  6  from  R  -  r  to  R  and  using  the  boundary 
condition  c^(R  -  t  )  =  o,  the  constant  N  is 


4nDeCi(R) 

[1/(R  -  T  )  -  l/R] 


Integrating  equation  6  from  R  -  r  to  r,  and  using  the  expression 
for  N,  the  interior  concentration  field  is 


ci (r)  = 


Ci(R) [ 1/(R  -  r  )  -  1/r] 
[1/(R  -T)  -  1/R] 


We  use  mass  continuity  at  the  droplet  surface,  i.e., 


4nR2De  a  c^/ar|  R  =  4irR2Dgac  eX/arlR 


where  Dg  is  the  gas  phase  diffusion  coefficient  of  ammonia  to 
determine  the  surface  ammonia  gas  concentration  c^(R).  With 
v2cex  =  0  and  cex  =  cex(  °°  )  as  r  —  °°  ,  the  exterior  solution 
becomes 


cex(r)  «  [cex  (R)  -  cev  (~)]R/r  +  cex  («) 


Substituting  equations  7  and  9  into  equation  8,  the  interior 
solution  is  written  as 


Cj^r)  = 


[R/(R  -  r  )  -  R/r]DgCex(  -  ) 
{De  +  [R/(R  -  r )  -  l]Dg j 


VlWivIviv! 


Porous  Shell  Diffusion-Controlled  Reaction  Model 


The  transport  of  ammonia  gas  through  the  porous  shell 
is  determined  by  substituting  equation  10  into  equation  8.  The 
molar  transport  rate  of  ammonia  gas  is 


4*RDeDgcex(  00  > 


Using  the  definition  for  £,  and  because  the  transport  rate  of 
ammonia  gas  is  equal  to  the  rate  of  production  of  ammonium 
phosphate,  the  porous  shell  diffusion-controlled  reaction  rate 
becomes 


d£/dt  = 


4ffRDeDgCex(  °°  ) 


np(°){De  +  [R/(R  ~  T  )  ~  l]Dg\ 


Using  equation  5  to  express  explicitly  t  in  terms  of  £ ,  the 
porous  shell  diffusion  controlled  reaction  model  for  a  layered 
acid  droplet  with  a  shrinking-core  is 


(De-Dg)(^-J,c)-2ncS  DgR{  l(R-tc)3-3np(0)(^-^c)/(47rc*  )  j2/3-(R-tc)2}/np(0) 


=  47rDeDgRCex(ac)t/np(0) 


where  £_  is  the  extent  of  reaction  at  crystallization. 


5.  COMPARISON  BETWEEN  EXPERIMENT  AND  THEORY 

Table  1  shows  the  list  of  conditions  (i.e.,  droplet 
size,  ammonia  gas  partial  pressure  time  of  crystallization 

tc,  and  extent  of  reaction  at  crystallization)  for  which  theory 
and  experiment  are  compared.  The  shell  thickness  at  crystal¬ 
lization  as  calculated  from  equation  3  is  also  shown  in  Table  1. 
For  these  calculations,  the  shell  porosity  e  is  assumed  to  be 
0.5.  The  shell  thickness  increases  with  increasing  droplet  size 
and  decreasing  ammonia  gas  partial  pressure.  For  rapid 


saw 
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reactions,  or  alternatively  small  droplets  and  large  ammonia  gas 
partial  pressures,  thin  shells  form,  while  for  large  droplets  and 
small  ammonia  gas  pressures,  thick  shells  form.  These  calcu¬ 
lations  are  supported  by  observations  with  a  telemicroscope  that 
permits  droplet  sizing. 

The  best  fit  between  theory  and  experiment  is  obtained 
by  adjusting  the  value  of  the  product  De(l-f)  so  that  the 
difference  between  data  and  theory  is  a  minimum  (see  Appendix) . 
For  the  shrinking-core  model  to  be  valid,  it  is  necessary  that 
the  phosphate-saturated  annulus  remain  saturated  during  droplet 
reaction.  The  characteristic  time  for  the  internal  liquid  acid 
core  to  become  uniformly  mixed  is  given  by  R2/D1,  where  is  the 
liquid  phase  diffusion  coefficient  of  ammonium  phosphate.  Using 
Stoke's  law  to  estimate  D-p  the  characteristic  time  for  uniform 
mixing  of  the  acid  core  is  about  300  sec.  Thus,  for  the  reaction 
times  of  interest,  the  saturated  annulus  remains  saturated  and 
the  shrinking-core  model  is  valid.  As  soon  as  this  zone  is 
depleted,  the  shrinking-core  model  is  not  valid.  The  extent  of 
reaction  at  which  the  saturated  zone  is  depleted  is  given  by 

*D  "  *c  +  4ffCp  (R  -  T  c)3  -  6  3  ] / np ( 0 ) 

where  &D  is  the  extent  of  reaction  at  which  the  saturated  annulus 
is  depleted.  For  run  46,  the  extent  of  reaction  at  depletion  is 
0.55.  At  this  point,  the  porous  shell/shrinking-core  model  is  no 
longer  valid.  In  the  following  analyses,  we  compare  the  semi- 
empirical  model  to  the  reaction  data  up  to  the  time  that  the 
saturated  zone  is  depleted. 

Figure  4  shows  the  comparison  between  experiment  and 
theory  for  the  four  separate  runs  listed  in  Table  1.  The  four 
reactions  correspond  to  reactions  conducted  at  increasing  ammonia 
gas  partial  pressures.  We  adjust  the  product  De(l-e)  to  minimize 
the  deviation  between  experiment  and  theory.  The  product 
decreases  by  an  order  of  magnitude  as  the  ammonia  gas  partial 
pressure  increases  by  almost  an  order  of  magnitude.  Assuming  a 
constant  porosity  of  0.5,  the  shell  diffusion  coefficient 
decreases  from  0.025  to  0.001  cm2/sec  as  the  ammonia  gas  partial 
pressure  increases  from  115  to  1000  dyn/cm2.  While  it  can  alter 
the  range  of  the  diffusion  coefficient  variation,  uncertainty  in 
the  porosity  e  cannot  account  for  the  trend  that  with  increasing 
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ammonia  gas  pressure  the  shell  diffusion  coefficient  decreases 
substantially. 
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Table  1.  Experimental  Conditions  for  Reactions  Between 
Phosphoric  Acid  Droplets  and  Ammonia  Gas. 


Run 

Radius  (jjm) 

pAg(dyn/cro2) 

tc(sec) 

S  c 

n) 

46 

32 

115 

20.0 

0.36 

1.01 

42 

33 

400 

7.0 

0.31 

0.99 

50 

30 

800 

4.5 

0.30 

0.69 

48 

29 

1000 

3.5 

0.18 

0.62 

According  to  Wheeler's  parallel  pore  model,2  the  porous 
shell  diffusion  coefficient  is  related  to  the  single  pore 
diffusion  coefficient  Dp  by  the  relation 


De  " 


where  e  is  the  porosity  of  the  shell  defined  as  the  ratio  of  the 
pore  surface  area  to  the  shell  surface  area.  The  tortuosity  k 
describes  the  deviation  of  the  pore  from  an  ideal  cylinder.  The 
diffusion  coefficient  decreases  if,  the  individual  pore  diffusion 
coefficient  decreases,  the  porosity  of  the  shell  decreases  or  the 
tortuosity  increases.  In  this  study,  we  assume  that  the  pore 
tortuosity  is  equal  to  unity.  The  Appendix  shows  that  the  porous 
shell  diffusion  coefficient  is  inversely  proportional  to  (1-e). 
This  fact  arises  from  the  reciprocal  dependence  between  the  shell 
thickness  and  the  shell  diffusion  coefficient.  Explicitly 
showing  the  dependence  on  e ,  equation  14  becomes 


DP  = 


e(l  -  f  ) 


where  A  is  the  value  of  the  product  De(l-e)  that  satisfies 
equation  13.  Recall  that  for  each  run,  A  will  assume  a  different 
value.  Since  the  product  e(l  -  e  )  does  not  vary  significantly 
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(less  than  a  factor  of  two)  as  the  porosity  varies  from  0.2  to 
0.8,  we  use  a  mean  value  approach.  For  f  =  0.5,  the  pore 
diffusion  coefficient  is  Dp  =  4A. 

This  relation  permits  the  evaluation  of  the  individual 
pore  diffusion  coefficient  from  knowledge  of  the  parameter  A,  the 
evaluation  of  which  does  not  require  an  explicit  knowledge  of  the 
shell  porosity. 

If  the  capillary  radius  is  much  less  than  the  mean  free 
path  of  the  gas,  then  diffusion  through  the  capillary  is 
described  by  Knudsen  flow.  The  individual  pore  diffusion 
coefficient  for  a  cylindrical  capillary  is 

Dp  =  4rc(2KLT/rrM)  1/2/3  (16) 

where  rc  is  the  individual  pore  radius,  K  is  Boltzmann  constant, 

L  is  Avagadro's  number,  T  is  the  temperature,  and  M  is  the 
molecular  weight  of  the  diffusing  gas.  As  the  capillary  radius 
decreases,  the  diffusion  coefficient  decreases.  Table  2  shows 
the  pore  radius  calculated  by  substituting  the  experimentally 
determined  pore  diffusion  coefficient  (equation  15)  into  the 
Knudsen  flow  diffusion  coefficient  (equation  16) .  As  the  ammonia 
gas  partial  pressure  increases,  the  pore  radius  decreases.  The 
dependence  of  the  pore  radius  on  the  ammonia  gas  partial 
pressure,  or  equivalently  the  reaction  rate,  explains  the 
phenomenon  of  encapsulation  that  is  observed  at  ammonia  gas 
partial  pressures  exceeding  1000  dyn/cm2 .  At  these  pressures, 
the  pore  size  could  be  so  small  that  the  ammonia  gas  molecule  is 
too  large  to  penetrate  the  individual  pore.  That  is,  the  shell 
acts  like  a  molecular  sieve.  Indeed,  the  spherical  equivalent 
radius  of  an  ammonia  gas  molecule  is  1.54  Angstroms,  and  Table  2 
shows  evidence  that  the  pore  radius  is  approaching  the  physical 
radius  of  the  ammonia  gas  molecule. 

6 .  SUMMARY 

We  derive  a  semiempirical  reaction  model  that 
describes  the  porous  shell  diffusion-controlled  reaction  of 
layered  droplets.  Specifically,  we  present  a  model  for  the 
reaction  of  phosphate-layered  acid  droplets  immersed  in  ammonia 
gas.  The  model  is  parameterized  with  an  effective  shell 
diffusion  coefficient  to  describe  the  transport  of  ammonia  gas 
through  the  phosphate  shell.  The  further  depletion  of  phosphoric 


acid  reactant  inside  the  droplet  is  modeled  after  a  shrinking- 
core  model.  We  compared  the  reaction  model  to  data  obtained  from 
the  electrodynamic  balance.  We  found  that  the  shell  diffusion 
coefficient  depends  on  the  droplet  reaction  rate,  decreasing  more 
than  an  order  of  magnitude  as  the  ammonia  gas  pressure  increases 
by  an  order  of  magnitude.  Furthermore,  the  average  pore  radius 
of  the  shell  decreases  with  increasing  droplet  reaction  rate. 

This  accounts  for  the  observed  negative  correlation  between  the 
shell  diffusion  coefficient  and  ammonia  gas  pressure.  For  large 
ammonia  gas  partial  pressures,  the  shell  diffusion  coefficient 
decreases  to  zero,  and  the  droplet  is  encapsulated.  Future  work 
will  be  directed  toward  the  determination  of  the  porosity  of  the 
phosphate  shells,  and  thereby  a  more  accurate  determination  of 
the  shell  diffusion  coefficient. 

Table  2.  Pore  Radius  as  a  Function  of  Ammonia  Gas  Partial 
Pressure. 


rc(A) 


PA^(dyn/cm2) 
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THE  DEPENDENCE  OF  THE  SHELL  DIFFUSION 
COEFFICIENT  ON  THE  POROSITY 


The  dependence  of  the  shell  diffusion  coefficient  De  on 
the  porosity  t  of  the  shell  is  derived  from  the  following 
considerations.  Using  the  definition 


3c  J[S(r)  -  l]r2dr/p 


we  can  rewrite  equation  3  as 


=  R  -  [R3  -  A/ (1  -  e  ) ] 1/3 


Since  r  <<R,  A/[RJ(1  -  e  )]<<!,  and  the  shell  thickness  reduces  to 


rc  =  R 


-  R|l  -  A/ [3  R3(l 


3R^ (1  -  e  ) 


Thus,  as  the  shell  porosity  increases,  the  shell  thickness  at 
crystallization  increases.  To  determine  the  dependence  of  the 
shell  diffusion  coefficient  on  the  shell  porosity,  equation  A  is 
substituted  into  equation  13.  Preliminary  to  this  substitution, 
we  use  the  approximation 


(R  -  rc)3»3np(0)  (  S-  £c)/(4rrc^) 


to  simplify  equation  13.  Approximation  equation  B  is  clearly 
valid  during  the  early  stages  of  reaction  when  £  =  £c.  Using 
this  approximation  to  expand  the  expression  in  equation  13,  we 
obtain 


( De  -  Dg)  (&-&c )  +  {^-^c)DgR/lR  -  ic)  =  47rDeDgRCexH+/nP(0)  (< 


With  r  <<R,  equation  C  becomes 

C 


=  4ffDgRcex(  oc  )  t/np(0) 


(I 


Thus,  holding  all  pertinent  parameters  constant,  the  ratio  rc/D( 
must  remain  constant.  Using  equation  A,  the  dependence  of  the 
shell  diffusion  coefficient  on  the  porosity  is 


De  =  constant/ (1  -  e  ) 
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